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SYSTEM AND METHOD FOR MODELING BIOLOGICAL SYSTEMS 

CROSS REFERENCE TO RELATED APPLICATION 
This application claims the benefit of piioiity of provisional U.S. patent application Ser. 
No . 60/304,289, filed July 9, 2001, which is incoiporated heiein by reference , 
5 BACKGROXJND OF THE INVENTION 

1 Field of the Invention 

The present invention relates geneially to modeling biological systems. More 
specifically, the present invention relates to a system and method foi" modeling biological 
systems using a component-based architectuie; the present invention also relates to a distributed 
1 0 computing aichitectuie for- implementing a system and method foi modeling biological systems . 
2. Description of the Background Art 

Software models of biological systems are extremely valuable to researchers. Such 
models allow I'esearchers, based on current theories about Uie innei workings of such systems, to 
predict the functioning of a biological system under- specific circumstances , Before the 
15 development of software models, die evaluation of detailed models of most biological systems 
was impractical, due both to the large number* of calculations that needed to be performed and to 
the need to undei stand thoroughly every aspect of the model . 

One example of a software model of a biological system is the computational model foi^ 
simulating the electrical and chemical dynamics of the heart that is desciibed in U.,S Patent No.. 
20 5,947,899 (Computational System and Method for* Modeling the Heart), which is incorporated 
herein by reference An example of a cell-level model of a biological system is the Virtual Cell, 
a software package developed at the University of Connecticut. The Virtual Cell and its 
capabilities is described in some detail in flie following references: JC. Schaff, B.M, 
Slepchenko, & L.M. Loew, "Physiological Modeling with the Virtual Cell Framework," in 
25 Methods in Enzymology, vol. 321, pp,. 1«23 (M. Johnson & L. Brand, eds.. Academic Press, 
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2000); J. Schaff & LJM. Loew, "The Viitual Cell," Pacific Symposium on Biocomputing 4: 228- 
39 (1999); C C. Fink et al,,, "An Image Based Model of Calcixmi Waves in DifFeientiated 
NeuioWastoma CeUs/' Biophys., J.. 79: 16.3-183 (2000); and J. Schaff et at, "A Geneial 
Computational Fiamewoik foi Modeling Celluiai Stractui^ and Function," Biophys J. 73(3): 
5 1135^46 (1997) The Viitual Cell software package is also desciibed in U.S., Patent No. 
6,219,440 (Method and Apparatus foi Modeling Cellular Stiuctuie and Fimction), which is 
incoiporated herein by i efeience 

Further examples of biological simulation software for modeling of biological and 
physiological systems inclnde: DBsolve (see I. Ooiyanin et al, "Mathematical Simulation and 

10 Analysis of Cellular' Metabolism and Regulation/' Bioinfoimatics, vol 15, pp., 749-58 (1999)); 
GBPASI (see P. Mendes & D. Kell, "Non-Lmear Optimization Of Biochemical Pathways: 
Applications to Metabolic Engineering and Parameter Estimation,^' Bioinformatics, vol, 14, pp, 
869-83 (1998); P. Mendes, "Biochemistry By Numbers; Simulation of Biochemical Pathways 
with GEPASI 3," Trends Biochem. Sci., vol 22, pp 361-63 (1997); P Mendes & D. B Kell, 

15 "On the Analysis of the Inverse Problem of Metabolic Pathways Using Artificial Neural 
Networks," Biosystems, vol, 38, pp, 15-28 (1996); P. Mendes, "GEPASI: A Software Package 
for Modeling the Dynamics, Steady States and Control of Biochemical and Other Systems," 
Comput. Appl. Biosci., vol. 9, pp. 563-71 (1993)); NEURON (see M. Hines, "NEURON: A 
Program for" Snnulation of Nerve Equations/* Neuial Systems: Analysis and Modeling (F, 

20 Eeokman, ed., Kluwer Academic Publishers, 1993)); GENESIS (see IM. Bower & D Beeman, 
The Book of GENESIS: Exploring Realistic Neural Models with the General Neural Simulation 
System, (2d ed., Springer^Veilag, New York, 1998)) . 

Numeions other simulation packages have been applied to modeling biological and 
physiological systems including: MetaCon (a DOS-based metabolic contrcl analysis program 

25 available at ftp://bmshuxleybrookes.ac ulc/pub/softwaie/ibmpc/metacon); Talis (a visual and 
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inteiactive real-time tool foi simulating metabolic pathways, gene circuits and signal 
tiansduction pathways); NetWoik (a Java applet foi intei active simulation of genetic netwoiks); 
SCAMP (a command-line driven software package lunning on the Atari ST and MS-DOS 
operating systems; capable of simulating steady-state and transient behavioi of metabolic 
5 pathways and calculation of all metabolic control analysis coefficients); MIST (a biological 
pathway simulation package innning on MS Windows 3,1); MetaModel (J^S -DOS -based 
softwaj^ package foi steady-state simulation of ^ metabolic pathways); SCoP (a commercial 
simulation program that can be used to simulate metabolic systems); CONTROL (a DOS-based 
software package that uses the Reder matiix method to calculate control coefficients fiom 

10 elasticity values); FIuxMap (a simulation package that calculates metabolic fhixes based on 
metabolite balancing); BioTbeimo (a simulation package that calculates the feasibility of 
individual pathway reactions based upon Gibbs free energy values and metabolite 
concentrations); BioNet (a metabolic flux analysis package); and the Matlab Simuliok and 
Statef low simulation packages. 

15 Traditional software models of biological systems typically compiised a system of 

differential and/or algebraic equations, which often had to be tianslated into numerical 
appioximations of those equations and then linked to a solvei to solve the system of equations 
and to graphics libiaiies for displaying the results of the simulation,. Early iinplementations of 
biological simulation software treated the entii© model as a monolithic system of differential and 

20 algebraic equations Hence, to make any change in a model, the software code typically had to 
be substantially rewritten. Given the lapid development of biological knowledge, such 
lepiogramming has pioven to be inconvenient and inefficient,. In addition, software designed to 
model a specific biological system could not be readily adapted to model a different but related 
biological system.. Therefore, a system using reusable software components, each representing 

25 some aspect of a biological system, would be desirable. 
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One appioach to solving this pioblem consisted of carviag up a large, complex model 
into a seiies of subparts ox sub-models, and then otganizing these sub models in a hieiaichicai 
stnictuie . Following this modulai' piogiaimning approach, a progiam would no longei consist of 
only one single part but lathei would be divided into seveial smaller paits which interact through 
5 procedui e calls and which form the whole piogiam» The main piogiam would coordinate calls to 
pioceduies in separate modules and woiild hand over appropriate data aspaiameteis, 

A modulai programming approach does not, however, completely solve the above- 
desciibed pioblemS:, and, in fact, creates some new ones. Foi example, many sub-models use 
global variables (whose values maybe changed by other sub- models), making it diEQcult, i£ not 

10 impossible, to understand how a particular sub-model would function or beihave without 
referencing the entire model. Indeed, for complex hieiarchical models, it may be difficult to 
determine whether a particulai* variable is used oi' referenced before its value is first defined. 
Moreover, the actual simulation results would often depend upon the order in which the sub- 
models were parsed.. Hence, sub-models caimot be neadily reused as components in anothei 

15 model without redesign, and modifications to a model still would require a deep understanding 
of the algorithmic details of the software. 

One approach to developing easily reusable software components is to adopt an object-- 
oriented arcMtectui^.. A traditional object- oriented architecture, however, usually uses 
inheritance to deJBne the relationship between components and sub-components, which is 

20 frequently inappr opriate in the context of biological modeling,. Inheritance is a mechanism that 
allows one object or' class (refened to as the child class or subclass) to be defmed by reference to 
another' object or- class (referred to as the parent class oi- superclass) That is, the child class is 
able to access the attributes and methods of the parent class without the need to explicitly 
redefine those attributes and methods for- the child class. 
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Oftentimes, in biological model, components and sub-components do not have a classic 
parent-child relationship, Foi example, a highei-level "patent" component such as a 
computational component foi* modeling a cell may comprise oi encompass lowei-Ievel **child" 
components such as a nucleus component, a cell-membiane component and an endoplasmic 

5 leticulum component Howevei, these subceliulat components may not shaie, and probably 
would not shaie, all of the attiibutes of the paient component Hence, it would be inappropriate 
to define the subcomponents by inheiitence . 

Anothei' diawback to conventional approaches to biological simulation is that they 
typically use a "single computei" aichitecture,. That is, that piocesses and calculations are 

10 peifoimed by a single pxocessox on a single machine. In some cases, models of complex 
biological systems can be so laige that they cannot be effectively evaluated on a single computei 
using piesently available haidwaie (oi can only be evaluated usiag prohibitively expensive 
haidwme). It would thciefore be useful to peifoim calculations lelating to such models in 
par allel on multiple computet^ in a distiibuted computing environment.. Hence, it would be 

15 desnable to ci"eate reusable software components capable of being used in a distributed 
computing enviEonment Moreover, because available computers are not always compatible, it 
would be desirable to utilize a platfoim- independent aichitectuie to allow the use of the 
piocessing capabilities of many types of computeis., 

SUMMARY OF THE INVENTION 

20 In accoidance with a fii-st embodiment of the invention, theie is piovided a system and 

method foi' modeling a biological system, comprising a pluiality of software components, each 
of which is connected to at least one other softwaie components, and each of which directly 
communicate only with parent ot child software components In accordance with anothei aspect 
of the invention, fheie is provided a system and method foi modeling a biological system, 

25 comprising a usei interface and a simulation engine, wherein tLe user inteiface resides on a 
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difierent computei than the simulation engiae . In vaiious embodiments, the simulation engine 
may include one or more of the following conrponents: an ODE solvei, a Solver Factoiy> and a 
Model Pactoiy,. In vaiioxis embodiments, the usei interface may include one or more of the 
following components: an XML paising module, a mathematical equation generation module, 
5 and a visualization/display module. 

In accoidanca with anothei aspect of the invention, there aie provided a method and 
system of modeling a biological system, wherein each of a plurality of instantiated software 
components is connected to at least one othei' of the plmahty of software components, each 
software component representing at least one of a physical poxtion of the modeled biological 

10 system and a functional aspect of the modeled biological system and data being directly 
conmmnic^ed only between software components directly connected to each other. Also 
provided is computei -readable mediimi having stored thereon computei- executable instructions 
for performing the above-described methods . 

In accordance with anothei* aspect of liie invention, there ai« provided a method and 

1 5 system for modeling a biological system, including a processor and a memory in communication 
with the processor, the processor- causing a pluiality of software components to be instantiated in 
the memory, the processor causing each of the plurality of software components to be connected 
to at least one othei- of the plurality of software components, each softwaie component 
r epresenting at least one of a physical portion of the modeled biological system and a ftmctional 

20 aspect of the modeled biological system^ and data being diiectly communicated only between 
software components directly connected to each other . 

In accordance wi"fli anothei aspect of the invention, there are provided a method and 
system of modeling a biological system, wierein a software component is defined and 
instantiated dynamically and the software component is subsequently redefined.. In one 

25 embodiment, the software component is defined by prototype . 
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In accoixJance with anoUiei aspect of the invention, thoi e is piovided a computet -loadable 
medium haviag stoied thereon a soflwaie component lelating to a model of a biological system, 
including an attrfbute relating to at least one state vailable, an attiibute lelating to an 
initialization method, and an attiibute lelating to a method of deteraiining the rate of change ovei 
5 time of the at least one state vaxiable. 

In accoidance with anothei aspect of the invention, thet e is provided a conq)uter-readable 
medium having stored thereon an adapter relating to a model of a biological system, including an 
athibute relating to an expected output variable of a first component, an attribute relating to an 
expected input vaiiable of a second component, and an attiibute lelating to a function foi 

1 0 mapping the value of the output variable to a value fliat is a valid value for the input vatiable., 

In accorxiance with another aspect of 1he invention, there are provided a method and 
system of evaluating at least one chaiacteiistic of a modeled biological system, including a 
plurality of connected softwax'e components, wherein the amount of time necessary to evaluate at 
least one chatacteiistio on a Sxst computei is estimated, a plurality of unselected connected 

1 5 software components is selected that can be evaluated within &e predetermined length of time on 
another computer" until the estimated time does not exceed a predetermined length of time, data 
relating to the attributes of any selected software components, data indicative of the state of the 
selected software components, data indicative of the rate of change of the state of the selected 
software components, and a direction to evaluate the selected software components without 

20 reference to the other software components in the modeled biological system are transmitted to 
the other- computer, any unselected software components are evaluated without reference to the 
other software components m the modeled biological system, data indicative of a rate of change 
are received from each other computer to which data were transmitted, and the entite modeled 
biological system is evaluated, 
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In accordance with another aspect of the invention, theie are provided a method and 
sj^em foi evaluating at least one chaiacteristic of a modeled biological system, including a 
plmality of connected softwaii& components, the system including a first processoi, a memoiy in 
coimnxmication with the Gist processoi, and at least one additional piocessoi in coimnunication 
5 with the first processox, the first processor estimating the amount of time necessary to evaluate 
the at least one characteiistic on a fiist computer; the first processor selecting a plurality of 
comiected software components that can be evaluated within the piedetennined lengfli of time on 
one of said at least one additional piocessor and 1±iat have not previously been selected if the 
estimated time exceeds a predetermined length of time, the first processor causing to be 

10 tr ansmitted to liie one of said at least one additional processor data relating to Hie attributes of the 
selected softwaie components, data indicative of the state of the selected software components, 
data indicative of the rate of change of the state of the selected software components, and a 
direction to evaluate the selected software components without leference to the other software 
components hi the modeled biological system if' any softwaiTe components have been selected, 

15 the first processoi evaluating any unselected software components without lefeience to the other' 
softwai'e components in the modeled biological system, the first processoi^ receiving fi:om the one 
of said at least one additional piocessor, data indicative of a late of change and the first processor 
evaluating the entne modeled biological system based on the received data and the calculated 
data. 

20 BRIEg DESCRIPTION OF THE DRAWINGS 

FIGURES lA, IB and IC illustrate systems in accordance with several embodiments of 

the present invention.. 

PIGURE 2A illustiates a non-hierai"chical view of the connections between multiple 
software components in a model in accoi dance with two embodiments of the piesent invention.. 
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FIGURE 2B illustrates a hieiarcliical view of the connections between multiple softwaie 
components in a model in accoi daace with two embodiments of the piesent invention. 

FIGURE 3 illustiates a method in accordance with a &st embodiinent of the pi-esent 
invention 

5 FIGURE 4 illustiates a method in accordance with the second emboditnent of the piesent 

invention. 

FIGURES 5 A, 5B and 5C depict several appioaches to implementing a distiibuted- 
component aichitectuie , 

DETAILED BESCRIPTION OF THE TNWTSmON 
10 The following teims shall have, foi the pmposes of this application, the lespective 

meanings set foith below, 

Adaptei: A softwaie component capable of adapting the inteiface of a software 
component to conform to the lequiiements of another' software component. An adaptei can 
change not only the iiqjuts and outputs of a software component, but its function as well, 
15 Examples of uses of an adaptei include: changiag the physical units of an output of a model 
component; changiag the mathematical problem type of a component; and adding equations or- 
inequalities to the systems of equations defined by a model component,. See Dan Becker, 
"Design Networked Applications iu RMI Using the Adaptei Design Pattern," JavaWoild (May 
1999), 

20 Black Box Component: A software component that hides all of the details of its 

implementation such that the only details that need be known about the component to use it are 
details lelatrag to its input variables and output variables, aud the publicly available functions the 
component peifoims . 

Chamcteiistio: A chaiacteristic of a modeled biological system includes, but is not 
25 limited to, any state variable of any component in the system A characteristic can also relate to 
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an, expeiimental paiametei' that is not embodied in any specific state vaiiable,, The evaluation of 
a characteiistic includes, but is not limited to, the determination of the value of a state vaiiable of 
any component in the system at a specific time, or aftei the occunence of an event (such as the 
addition of a chemical to the system), oi' the determination of the late of change of the value of 
5 such a state variable* 

Definition by Prototype: The definition of a software con^nent oi' othet software object 
by the description of a particular' instance of that software component or other' software object, 
r ather than by means of an abstract declaration, 

Dynamic Definition: The definition of a software component or other softwaie object oi 
10 other variable in such a fashion that any later redefuiition of the software compojnent or other 
software object or other' variable will cause any previously instantiated softwaie component or' 
other software object or declared variable to conform to the revised definition. For example, a 
dynamically defined cell membrane component can include a method that defines the rate of 
sodium transport across it. That method can be used to calculate ttansport in the early stages of a 
15 simulation, but be redefined by other code to calculate transport differently. The redefined 
method would be used to calculate transport after such redetinition., 

Fimctional Portion: Any portion of a modeled system that is not a physical portion of 
that system, such as a process for" absorption of a chemical through a cell membrane. 

Input Variable; A variable stored in a software component as a property or other field 
20 theoceof that can oi" must be passed to that software component at initialization or- later,. 

Instantiation: The creation of a specific instance of a software component or other' 
software object at runtime. 

Output Variable; A variable stored iq a software component as a property or other field 
thereof that can or must be passed fiom that software component to another- softwar e component 
25 01 other object or fdnction 
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Physical Poition; Any portion of a modeled system that is a tangible sabpait of liie enth"© 
system, such as a cell, a cell membiane, or a cell nucleus . 

Software Component: Within the context of the present invention^ a software component 
is a xieusable progiam building block that can be combined with othei' components in the same oi 
5 other computers in a distiibuted network to foim aa application, Examples of components 
include: a single button in a ^aphical usex' inteiface, a small inteiBst calculatoi; an intei&ce to a 
database managei. Components can be deployed on diffeient seivers in a network and 
commimicate with each other for" needed services^ A component nms within a context called a 
contaiaer. Examples of containers include pages on a Web site, Web browsers, and word 

10 processors. Software components can also represent entities in the real-world, such as models of 
biological processes and sub-processes , They can also represent moie abstract components 
necessary for building, storing, simulating and visualizing models An example of an abstract 
component is an Adapter, which was defined above. 

State Variable: A property oi other field of a software component oi' other software 

15 object that represents a physical oi' functional aspect of a model, For example, a variable 
(whether oi" not publicly accessible to at least one othei^ object or method) representing the 
concentration of a particular chemical in a cell is a state variable. AU ii^ut and oulput variables 
are state variables, although a state variable need not be an input or- output variable . 

Commonly owned US, Patent Number 5,947,899 titled "Computational System and 

20 Method for* Modeling the Heart" and commonly owned U S . Patent Application Serial Numbers 
09/295,503 ("System and Method foi" Modehng Genetic Biochemical, Biophysical, and 
Anatomical Information: In Silico Cell"), 09/499,575 ("System and Method for Modeling 
Genetic Biochemical, Biophysical, and Anatomical Information: In Silico Cell"), 09/599,128 
("Computational System and Method for Modeling Protein Expiession")^ 09/723,410 ("System 

25 for Modeling Biological Pafhways")^ and 09/898,151 (Method and System for Modeling 
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Biological Systems) piovide fuithei discussion and explanation of vaiious aspects of the 
modeling of biological systems and are hereby incoipoiated by reference in theii entirety. 

Refening to Figure lA, a block diagram of a system in accordance with a fiist 
embodiment of the piesent invention is illustsated, Computer 100 includes processoi^ 102 and 
5 memory 104 connected to processoi 102, Computer 100 can be a maiofiame computer, 
minicompntei, rcdcrocomputei' (including a personal computet oi a woikstation), oi othei' 
computing device,. Ptocessoi 102 can be an Intel Pentium oi equivalent miciopiocessoi, 
although othei suitable processoi s can be utilized, Memory 104 can be permanent stor age, such 
as a hard drive, temporary memory, such as random access memory, or a combination of 

10 permanent storage and teniporaiy memory Memory 104 typically includes both one or more 

, hard drives and a form of random access memory The fust embodiment is most suitable for 
implementations in which a user 's computer is likely to be powerful enough to evaluate models 
likely to be used by a typical user within a reasonable amount of time At the pr esent time, the 
second and third, embodiment described below with respect to Figures IB and IC are more 

15 appropriate for the evaluation of complex models; howevei, if the processing capabilities of low 
cost computers continue to improve, the first embodiment could become more appropriate for the 
evaluation of complex models., 

Refening to Figure IB, a block diagram of a system in accordance with a second 
embodiment of the piesent invention is illustrated Computer 100a includes processor 102a and 

20 memory 104a connected to processor 102a,. Computer- 100a, which functions as a server, can be 
a mainframe computer*, minicomputer, microcomputer (including a personal computer' or a 
workstation), or other computing device, but is typically a supeicomputei, such as SGI machine, 
running the Oiigin 2000 operating system, Processor 102a can be a RIOK or' equivalent 
processor-, although other suitable processors can be utilized. Memory 104a can be permanent 

25 storage, such as a hard drive, tempoiaiy memory, such as random access memory, or- a 
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combimtion of peimanent storage and ten^oiaiy memoiy. Memory 104a typically includes 
both one oi moi e bar d drives and a foim of landom access jnemoiy . 

Computei 100b includes processoi 102b and memoiy 104b connected to processoi 102b . 
Computei 100b, which functions as a client, can be a mainfiame computei; minicomputei, 
5 miciocomputei (including a peisonal computei of a woikstation), oi othei computing device, but 
is typically a personal computer, such as an IBM compatible personal computei lunning the 
Miciosofi: Windows NT opeiating system., Piocessoi' 102b can be an Intel Pentium family oi 
equivalent microprocessoi, although othei suitable processors can be utiEzed, Memoiy 104b can 
be permanent storage, such as a haid diive, temporaiy memoiy, such as random access memoiy, 

10 or a combination of permanent storage and temporary memory* Memoiy 104b typically includes 
both one oi moie hard drives and a form of random access memoiy . Computei 100a can 
communicate with computer 100b across the Internet, a private network, by direct modem to 
modem communication across telephone lines, or by other methods. 

Referring to Figure IC, in a third embodiment of the present invention, a plurality of 

15 computers 100a through lOOn are utilized and connected together by the Internet oi^ other 
network.. Large models can be broken down into multiple portions, with each portion being 
evaluated on a separate computei, with data from the evaluated portions being recombraed and 
analyzed on a single solvoi- computer: In yet other embodtmeats, this analysis can itself be 
broken down into portions and each portion can be analyzed on a separate computer^ with a 

20 master solver' computer combining the analysis to reach a final result. The use of multiple 
computers in connection with biological modeling is described in fuithei detail in connection 
with the method illustiated in Figure 4 below. 

Refeiiing to Figuies lA through IC, biological model softwaie 106 is stored in memoiy 
1 04 In the second embodiment, part of biological model softwaie 106 is stored in memoiy 1 04a 

25 and part is stored in memoiy 104b. For- example, the user interface portions of 1iie softwaie, 



wo 03/007217 PCT/US01/S0491 

... 14 .. 

including the graphics fimctions, can be stored in memory 104b at the cJient, while the portions 
of the software dealing with the evaluation of models can be stored in memory 104a at fee 
server. Biological model software 106 can be written in C++, Java, or one or moi« of many 
alternative programming languages utilizing object oriented methodologies. Numerous 
5 references, including Erich Gamma et al.. Design Patterns: Elements of Reusable nhf «..t- 
Oi^goted Software (Addison- Wesley 1995) and Bertiand Meyer, Obiect-Qriented Sofhyam 
CQ;nstmctiott (2d ed, Prentice-Hall 1997), teach how to program using an object oriented 
language. 

In a first, second, and third embodiments of the present invention, biological model 
10 software 106 preferably lepieseiits pojtions of a model (such as a single cefl, a gioup of cells, a 
cell nucleus, or a paiticulai pathway) as software components, In the third embodiment of the 
preset invention, the con^Kinents are prefeiably constiucted in accordance with a platform and 
language independent standard, such as die Common Object Request Brokei' Architecture 
(heremaftei "CORBA"), to allow the components to be transmitted to a wide range of other 
15 computers, such as computeis lunniag mcompatible operating systems. The CORBA 
specification is available fiom the Object Management Group of Framingham, Massachusetts 
(www omg org) and is hsidby mcorporated herein by reference in its entiiety . 

Model data 108 can be stored m a database, such as a relational or object oriented 
database, a flat file, such as an XML file, or other type of file In certain embodiments, data 
20 relatmg to adapters can be stored in the same location or in a separate database or file Oracle 8i 
or any suitable custom or commercial off-the-shelf relational, object oriented, or oflier database 
can be used to store model or' adapter data A portion or all of biological model software 106, 
model data 108, or both biological model software 106 and model data 108 can optionally be 
stored on an internal medium, such as a hard drive orxandom access memory chip or module, or 
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on a removable medium, such as a floppy disk, tape, CD-ROM, DVD-ROM, removable haid 
drive, oi other form of magnetic, optical, or other storage,, 

A variety of modeling languages can be utilized in accordance the present invention. 
One such language is CellML, an extensible marlaqj language designed for biological modeling 
5 Portions of the biological simulation software may also be coded using MathMI^ which was 
originally developed as hypertext markup language for displaying mathematical equations but 
has been adapted &i use in genemting code for numerical oonq)utatioii (rather than merely 
rendeiing the equations foi' screen display or printing) . 

Refening to Figures 2A and 2B, two views of the connsctions between seveial 
10 components are illustrated The non-hiesraichical view of Figme 2A illustmtes the comiections of 
component 200 with its sub-components, Figure 2A does not iUustiate all of the connections 
between the illustrated component and sub- component, and other components and sub- 
components, and the iHustrated component and sub-components ai« fer less complex than typical 
modeled con^jonents Component 200 has one connected input variable, T. which is provided 
15 by a component that is not shown, and erne connected output component, E, which is an output 
variable of child component 220, Component 200 provides six input variables to child 
component 220. namely its input variable T, four of its local variables, vj, vi, Cas, and Ca,,, and 
an output variable of child component 250, ib. Na and receives one output variable, E, fiom 
component 220, and one output variable, ib. ca, fiom grandchild component 230 Grandchild 
20 component 230 has two mput variables, E and Eca, both of which are local variables generated 
by child component 220, and one output variable, ib, which it provides to child component 
220. Grandchild component 240 has no connections, while child conqjonent 250 has no 
connected input variables and only one connected output variable, i^, na, which is a variable 
generated locally within child component 250 Figme 2B illustrates the same component, sub- 
25 compon«rt3, and connections in a hierarchical view 
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Not all input and output variables need be connected in a model , For example, input 
vatiable 202 and output vaiiable 204 of component 200 are not connected. In this case, a default 
initial value of iiq}ut vaiiable 202 is geneiated (if it is necessary) by some initialization procedure 
by component 200. 

5 As will be seen firom Figures 2A and 2B, in the exemplary embodiment, aH 

communication is conducted between comiKsnents that are diiectly connected to each othex . 
Limiting communication to diiectly connected components eliminates the need of components to 
know details relating to the internal functioning of other components. Even if a variable is 
passed from one component to a second component^i which may not use that vaiiable internally, 

10 and on to a tbiid component, the first component need not know what component will use that 
variable's value and the Ihird component need not know the somce of the value of that variable. 
Thus, components can be easily reused in a completely transparent jEashion. For example, a 
researcher desiring to evaluate a system consisting of five components representing cells A 
through E can take a model that consists of cells A through D and cell F, disconnect and r emove 

15 the component representing cell F, and add and connect a component foi cell E. Because all of 
the ioteiaotions between the components lepiesenting the cells aie handled on a direct 
component to component level, none of tiie coxmections between cells A through D among 
themselves aie affected Moreover', the creation of components composed of other" components 
is gieatly simplified Finally, these components can be connected without any knowledge of the 

20 detailed structur e of these components. 

Notably, a group of connected components can be viewed as a component itself ' In fact, 
a preferred method foi' creating a new component by assembling pr eviously created components 
is a technique known as composition. Like class inheritance, compositioa is a mechanism for 
adding new functionality to a software class In software design, composition is cunently 
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favored in most cases ovei inheritance because composition is moiQ flexible and does not expose 
the internal implementation of a class . 

The Composite Design Pattern, see Eiich Gamma et al., Design Patterns: Elements of 
Reusable Object-Oriented Software (Addison-Wesley 1995), is also useful fox assembling 
5 biological models fi om simpler models. It is veiy convenient to consider a model of biological 
processes to be made up of simpler processes., Foi' example, the behavioi of a livei cell can be 
conveniently described by the stmctinra of the cell and fhe pxocesses of the cell. In turn, it is 
convenient to considei even moie complex models, e g. a livei oigan to be made up fiom a 
hierarchy of complex models, e.g. liver cells. The Composite Design Pattern allows fhe 

10 software to automatically tteat any simple or complex model as a model component. 

Refeiiiag to Figme 3, a method in accoidance with an exemplary embodiment of flie 
piesent invention is illustrated In step 300, a plmalily of softwaie components^ such as those 
described above in conneGtion with Figures 2A and 2B, aie instantiated. The consonants can be 
eithei- (i) first declaied and subsequently instantiated (e.g., fhe component can be defined in a 

1 5 declamtions section and instantiated in fhe body of the piogram) oi (it) defined by prototype, that 
is to say (with respect to the first instance of each class of component) simultaneously declaied 
and instantiated (i.e , the component is both defined and instantiated at fhe same point in the 
pmgiam code),. Each component preferably comprises one or more (but typically many) state 
variables, which describe important chaiacteiistics of the components. A state variable can be an 

20 input variable, an output variable, both an input and an output variable, oi a purely local variable . 
The component pr efeiably comprises one or- mor e initialization methods for ioitiaHzuig any state 
variables that are not iaput variables,. The component preferably also comprises one oi more 
methods for' determining fhe rate of change of its state variables over* time 

In step 302, with respect to each coimection between an output variable of a first 

25 component and an iaput variable of a second component that is to be made in a model, it is 
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optionally determined whether an adapter is uecessaiy An adapter is necessaiy if the format of 
the output of the fiust component does not match the format of the input of the seicond 
component The use of muitq>le adapters can be necessary wiUi respect to some components, 
especially components having many input and output variables. For example, if the first 
5 component generates an output in degrees Celsius and the second component lequiies an input in 
degrees Kelvin, an adq)ter' is necessary to translate fiom degrees Celsius to degrees Kelvin,. 
Similarly, an adapter is necessary if data generated by the fir^t cornponent is slightly different 
fiom that requiced by the second component For example, the first component could gene:ate 
an output in units of for ce (such as pounds) while the second component lequired an input in 

10 units of mass (such as fcQograxns). An adapter would be necessary to map the first value to a 
second value using a formula including constants and a variable (the gravitational force of the 
eaith. or other relevant object), In some embodiments of the piesent invention, adapters ai© not 
utilized and the value of each output variable must be a valid value for the input of the input 
variable to which it is cormected., 

15 In step 304, an appropriate adapter is optionally automatically located and retrieved in 

each case in which an adapter is necessary . Ihe use of consistent nomenolatiire or stiong typing 
greatly facilitates the automatic retrieval of adapfters For example, if the output variable of a 
first component is of a degrees Celsius type or labeled as degrees Celsius, and the input variable 
of a second component is of a degrees Kelvin type oi labeled as degrees Kelvin, it is a simple 

20 matter to query a database indexed on two type or label fields (one relating to the input to the 
adapter and one relating to the output) to locate an appropriate adapter,. Alternatively, the user' 
can be required to identify an appropriate adapter- in embodiments in which adapters are utilized. 
In addition, in embodiments in which adapters are utilized, the user^ can be given the option to 
ovenide the automatic selection of an adapter- with another*, more appropriate adapter. 
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In step 306, if adapters aie used^ they are connected to tho appiopiiate components. In 
step 308, the components (oi adapteis connected to components) aie connected togethei by 
associating input vaiiables of components with output variables of ofliei compon^ts . 

In step 3 10, a leal ordering of the oidei of initialization of Ihe components of a model is 
5 CI eated if the initialization methods of any state variables depend on values of state variables of 
othei components. If a real ordering is not cieated, the ordei in which the components aie 
initialized can affect the values of state variables that depend on values of state variables of other 
components* A i^al ordeiing can be created by first creating a diiected graph comprising all of 
the dependencies of state variables on state variables of other components at the time of 
10 initialization and then usizig the partial oider sort algorithm to convert the directed graph to a real 
ordering . In step 312, each component in the model is initialized, 

in step 314, in embodiments in which dynamic declaiation is utilized, components are 
optionally redefined after- instantiation of specific iastances.. For example, if the intt eduction of a 
particular chemical causes a change in the basic functionit^ of a component, the component 
15 might be redefined upon the introduction of that chemical, Each previously defined instance of 
that component, as well as subsequently defined instances of the component^ would assume all of 
the characteristics of the component as redefined. 

In step 316, a plurality of components is optionally combined into a single reusable 
component. For example, biological model software 106 can include a method for automatically 
20 generating a wrapper allowing the user to view the plurality of components as a single 
component while ti eating the plurality as separate connected components for intemal purposes., 

Referring to Figme 4, a method in accordance with a second embodiment of the present 
invention and usable in conjunction with the method described above in connection with Figure 
3 is illustiated. Prior to the performance of ^ step 400, a model of a biological system is 
25 constructed using a plur ality of connected components and the characteristics to be evaluated ai^ 



wo 03/00721 7 PCT/USOl/50491 

-20- 

prpvided to biological model softwaie 106* in step 400, the amount of time necessaiy to 
evaluate flie chatacteristic or chaiacteiistics of the modeled system is estimated. In certain 
embodiments of the present invention, the amount of time need not be measmed with pi ecision 
because the estimate is used only to deteimine whethei multiple computeis should be used to 
5 evaluate the model (or optionally the number of computeis that should be used)., The estimated 
time can be deteimined using seveial methods.. The estimate can be calculated using only the 
numbei of software components used in the model, using the number of software components 
and the length of simulated time that the model will rim, usiog the mmibei of^ software 
components and the oider of magnitude (eg., n, n log n, n^ etc.) of each algorithm in the model 

1 0 (oi only selected algoritbms), or using other methods . 

In step 402, if the estimated time is longei than a predeteimined length of time, a plurality 
of connected software components is selected, The selected components must be diiectly 
connected to each ottier, alfhough each selected component need not be directly connected to 
each other- coirq)onent The selected con[5)onents must also be capable of evaluation within the 

15 predeteimined lengfia of time on another computer, The other computer can be any type of 
computer fiom a personal computei to a super computer ; hence, it is desitable to know the type of 
computei that will be used before the performance of step 402. The components can be selected 
by dividing the model into an appropriate number of pieces and may take the capabilities of each 
other- computer to which portions of the model will be sent into account or may take into account 

20 only the capabilities of the least powerftil of such cotnputer. The dividing lines between the 
portions of the model need not be logical. 

In step 404^ if any components have been selected in step 402, data relating to such 
components is sent to another computer together* with a direction to evaluate such components. 
The data should include at least data relating to the curr ent collective state of the components and 

25 theii' cmient collective rate of change . A platform independent standard such as CORB A can be 
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used to send the data to the othei computei . The collective state of the selected components can 
be l epi esented as the union of all of theii state vaiiables and then collective tate of change can be 
l epresented as the late of change of that imion. 

In step 406, steps 400 through 402 are repeated until any nemaining softwaie components 
5 can bo evaluated within the predeteimined length of time bi step 408, the first computei 
evaluates any remaining software components without reference to those sent to other 
con^xiters,. In step 410, data is received fiom the computers to which data was sent in step 404., 
In an exemplary embodiment, data indicative at least of a rate of change is received fiom each 
computei to which components were sent in step 404. In step 412, tiie entire modeled system is 
10 evaluated based on the results of steps 408 and 410 , 

Optionally, step 410 can be performed on a plmality of computers m a similar manner in 
complex models where the evaluation will likely consume an excessive amount of time, In this 
option, functional domain decomposition is utilized to perform various tasks comprising the 
evaluation of a model in paraileL. Ihe use of algebraic domain decomposition is well known to 
15 those skilled in the art. For example. United States Patent Number 5,659,788, which patent is 
incoipor ated by reference herein in its entnety describes algebraic domain decomposition, which 
is typically utilized to divide a representation of a geometric ol:)ject into subparts and to process 
such subparts in parallel. Functional domain decomposition involves the use of similar 
techniques in dividing the computations required to calculate the result of a function into tasks 
20 and to process such tasks ia paiallel. 

First, the process of evaluation must be divided into a gioup of tasks that can be 
performed m parallel, Depending on the particular model, a greater- or lesser degiee of 
parallelism can be present,. In any event, the final calculations lequiied to generate an answer' 
will in most cases not be capable of performance prior to the generation of olher- data.. Many 
25 earlier calculations, however, can often be peifonned in parallel . Open MP or MPI can be used 
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in this step, Second, at least some of the parallel tasks aie sent to othei computeis to be 
peifoimed and the insults aie leceived fiom sxicli computer . Third, the results of the paialiel 
tasks are combined. Fouith, any tasks that aie not paiallel aie peribimed by the fiist computer 
The analysis is then complete,. 
5 In an exemplaiy embodiment, a separate software componeait referied to as a solvei is 

instantiated to peifoim a single task oi^ gioup of tasks necessary to complete the evaluation of a 
modeL In the exemplaiy embodiment, the several solveis instantiated to perform the tasks 
necessaiy to evaluate a model can be diffeient from each other ia many ways . They can be 
designed to apply diffeient algorithms, run on different compfuting platforms, etc,. The factoiy 

10 method (also known as a viitual constructoi) is used to define an interface for creating a solvei' 
component, but allow subclasses to decide the class of solvei- component to instantiate . 

The software components described above can be implemented using a vaiiety of data 
structuiesn One way in which components of models and adapteis can be inq>lemented using 
XML.. Ofhei non-XML data stEUctuiiesai^ also suitable. 

15 A Distributed Computing Approach 

Anothei advantage of a reusable software -component arcMtectme is that it lends itself to 
implementation of distcibuted computing solution. In a traditional design of a biological 
simulation model, aU software components (except possibly foi the database) would reside on a 
single machine, typically a desktop computer . The drawback to this design is that, as a model 

20 becomes more and moie complex, die computational requirements of peifoiming the simulation 
will gradually overwhelm a single computer's capacity in terms of computing power and 
memojy At some point, it will become necessaiy to move the software application to a moie 
powerful computer- with the ability to handle the increased computational loads 

In addition to the scalability and peifoimance issues, anothei drawback is the need to 

25 leinstall the enttie software packages if one component, such as the solver, is upgraded and needs 
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to be replaced, Henoe, thei© is an inherent maintenance problem lelating to deployment of 
upgrades., 

The solution to these pioblems is to adopt a distributed conaputing architectuie, with the 
MSGT inteiface (UI) components and the computationally intensive solver components lunning on 
5 sepaiate machines,, Foi example, the UI can imx on the usei*s desktop compntei' (oi poitable 
conoputei), while the simulation engiae runs on a more poweiful servei, such as a supercomputer 
01 a cluster' of Linux machines . Optionally, if the problem's solution can be parallelized, there 
can be multiple numetical engines ixmning on multiple computers, each solving subpioblems of 
the larger pioblem. The IJI components may also be split up to nm on seveial machines . The 
10 distiibuted computing environment may comprise a number of coniputers directly linked in an 
xntianet, ot remote computers may communicate ovei' the internet or by other means, 

Theie aie many standard distiibuted-computing technologies/aichitectuies, including 
CORBA, DCOM, Enterprise JavaBeans, Remote Method Invocation (RMI), EJB, Parallel 
Virtual Machine (PVM), Message Passing Interfece (MPI) and OpenMP* Preferably, the system 
15 is implemented using CORBA, which is platform-independent, language-independent, network- 
pi otocol- independent and object-oriented The CORBA specification therefore allows the 
creation of intetchaageable, reusable software components in heterogeneous computing 
environmeaats. In comparison, the other standards suffer from one or moie drawbacks. For 
example, DCOM, although in theory platfoim-independent, is not well-supported on non- 
20 Microsoft machines and is not an open standard; Enterprise TavaBeans, RMI and EJB are 
language-dependent; and PYM, MPI and OpenMP, which aie parallel computing technologies 
for low-level message passing in a distiibuted environment, cannot diiectly ttansrrut software 
objects and theiefoi'e cannot be easily adapted to implement an object-oriented system., 

Because C^ORBA's Interface Definition Language (IDL) is independent of the 
25 implementation language, it is possible to choose the most appropriate piogiamming language 
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and execution enviionment foi each componeat of the system mder development . Moieovei , 
because all objects interact thiough the DDL interface, legacy code (siich as pre-existing 
numerical solveis, for example) can easily be integrated into fhe application. Alternatively, 
CORBA's dynanuc interfaces (BE and DSI) can be used instead of IDL 
5 Furtheimoiej because all communication between objects and clients aie enabled oi 

mediated by the Ob|ject Request Broker (ORB), CORBA providea interoperability among 
components designed to run on different machines. That is, a client can invoke an object on a 
remote computer tianspatently and without regard to the location of the object or operating 
system under which the object is running Finally, anolhez advantage of COKBA is that there 

10 exist software tools for implementing felly seemed distributed systems and for designing load- 
balanced applications. 

Figure SA depicts one embodiment of a distributed computiag aichitectaie, wher ein the 
simulation engine or ODE Solver" 510 is implemented as a distributed component running on a 
different host machine than the user interface oi UI 501 The UI 501 will typically run on the 

15 user's local deslctop computer^ and the ODE Solver 510 will typically run on a remote server' 
connected to user's machine on a local area network (LAN). Preferably, the UI 501 is 
iiiq>lemented in JAVA. The ODB Solver 510 is any numerical solver capable of solving a 
system of ordinary differential equations (ODEs). The HI 501 allows fhe user' to create, modify 
or' view a biological model, as well as to enter- and retrieve data relating to the model 

20 Prefierably, the UI 501 comprises a graphical user interface (GUI), and includes an editor (more 
preferably a WYSIWYG editor) for entering and viewing the underlying mathematical equations 
for the biological model, as well as the other- modeling components or elements. For example, 
the UI 501 may include a "pathway editor*' that allows the user to create, edit and view pathway 
diagram representations of the biological system being modeled; oi' the UI 501 may incliide a 

25 "cell editor" that allows users to create, edit and view cell-level components for cell modehng. 
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In a prefen^d embodment, the simulation model is stored in an XML format, such as MafhML 
or CellML 

The tll 501, ra the embodiment depicted, haclttdes an XML Parsing component 502 and 
Math Equation Generation component 503 foi translating the user input into a set of equations in 
5 a foimat that can be solved by the ODE Solver 510. The UI 501 preferably includes a 
Visuahzation component 504 that allows the user to view the model and data, as well as the 
simulation results 

As depicted in Figuie 5A, the CORBA IDL solvei inteiface encapsulates the simulation 
engine (or ODE Solvei 510) and allows the UI 501 to access the ODE Solver 510 only fhiough 
10 the solver interface. In a prefeiied embodiment, the solver interface is implemented in C^-\- and 
is wrapped around a FORTRAN ODE numeiical solver. The remote server application activates 
a servant object that implements the solver interface; and the seivei ORB converts the object 
reference into a standard string format. The ORB on the UI side converts this string back into an 
object reference.. 

15 After the simulation model is constructed, the UI 501 acts as a cheat and uses the object 

reference to make a series of invocations of the remote solver object to simulate the model,. In 
response to each invocation, the solvei' computes a nurneiical solution and letmns the solution to 
the UI 501 for display . 

The embodiment depicted in Figure 5 A is perhaps the simplest design for separating the 

20 UI from the simulation engine, thereby uiiloading the computationally intensive components oft 
the user 's desktop machine . One drawback of this design is there is no means for the chent to 
instruct the remote server* to instantiate a new solver- object. 

A more flexible implementation is illustrated in Figure 5B, which is similai' to the 
embodiment depicted in 5A except that it introduces a Solver Factory object 520 The Factory 

25 Method (also referred to as a 'Viitual constructor'' or' a "creational design pattern") is a method 
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for absttacting the instantiation piocess. The Factor Method allows a class to defei instantiation 
to subclasses. That is, it defines an interface for creating an object but lets subclasses decide 
which class to instantiate . 

As shown in Figure 5B, the Sdlvei Factory 520 allows fo: the creation of a solver object 
5 on demand. The leinote server activates arid publishes the Solver Factory object. The server can 
then defer the instantiation of a solver* object until the client expHoifly requests it. To actually 
create a new solver object, the UI 501 invokes an operation on the Factory object, which then 
creates a new solver object and returns a reference for the new solver- object to the UI 501 . The 
client need not know where the actual solver object r esides; and the UI 501 needs only an object 

10 reference to initiate the simulation , 

Although the embodiment illustrated in Figure 5B is an improvement over the simple 
design show in Figure 5A, it is still iaefiicient. Notably, dming the numerical iategiation 
process, the ODE Solver 510 needs to know the right-hand side (RHS) ftmotions for the ODE 
system being solved. Because the RHS functions r eside on the UI machine, the ODE Solver' 510 

15 server must initiate a callback each time to retrieve the RHS functions., Since Ods callbaok event 
ocouis many times during Ihe simulation process, it may flood the commimication chaimel 
between the computation server and the UI machine, thereby degrading the perfoimance of the 
whole system,. 

Figure 5C illustrates an embodiment, improving upon the designs shown in Figures 5 A 
20 and 5B . In the embodiment shown in Figure 5C, the callback process in removed entirely, and 
the model generation process (Model Factory 530 and Model Generation Logic 540) is 
encapsulated into the IDL interface and implemented on the remote server,. In this embodiment, 
the UI 501 is a pure client and sends the model (pr eferably iu the form of an XML file) to the 
remote server through an IDL operation. The UI 501 invokes the solver' operations as before, but 
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only once for eaoh simulatioa. The simulation results are returned to Uie UI 501 for visualization 
only aftex the entiie simulation process is completed. 

As a i^ult, the commumcation-oveihead between the client machine and the servei are 
leduced significantly* Hence, this implementation is more efficient. Moieovei', multiple models 
5 may be simulated simultaneously using this implementation^ theieby making this implementation 
moie flexible than the two piioi implementations described in Figures 5A and 5B. 

All publications, patents, and patent applications mentioned in this specification Bxe 
herein incoipoiated by lefeience to the same extent as if each individual publication oi' patent 
application were specifically and individually designated as having been incoipoiated by 
10 reference. 

While this invention has been described with an emphasis upon piefeired embodiments, 
it will be obvious to those of oxdiuaiy skill in the att that variations in the pi efen'ed devices and 
methods may be used and that it is intended that the invention may be pxacticed otherwise than 
as specifically desciibed heiein Accoidingly, this invention iochides aU tnodifications 
15 encompassed within the spirit and scope of the invention as defined by the claims that follow.. 
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CLAIMS 

That which is clanned is: 

1 A method of modeling a biological system, compiising: 

(a) instantiating a plxiialily of software components; and 
5 (b) connecting each of the plumlity of softwaie components to at least one 

othei' of the plmality of software components, 

wheieia data are diiectiy commmiicated only between software components 

diiectly connected to each other, 

2 . The method of claim 1 , wher ein step (b) comprises associating one or moie irfpnt 
10 vaiiables of a software component with one or moie output vaiiables of anothei softwaie 
component. 

3 The method of claim 2, wherein the one oi" more iaput variables and the one oi 
mote output vaiiables each comprise at least one state vaiiable. 

4 . Ihe method of claim 1, ftuther compiising: 

15 (c) combining a plmality of connected components into a single reusable 

component 

5 . The method of claim 4, wherein the single i eusable component is a black box 
component 

6 The method of claim 1, ftuther" comprising; 
20 (c) initializiQg at least one field of each of flio plurality of connected 

components.. 

7. Ihe method of claim 6, wherein the vahie of the at least on© field in at least one 
of the plurality of connected components depends on the value of at least one field in at 

least one othei of the pluiality of connected components . 
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8, The method of claim 7, wherein with respect to each vahie that depends on the 
value of at least one field in at least one othei connected component, data relatiog to the 
dependency is passed to at least one othei' connected component . 

9, The method of claim 8, fetfhei' compiising: 

5 (d) cieating a leal oideiing representing the order in which the fields will be 

initialized 

10, The method of claim 2, further compiising: 

(c) connecting an adaptei' to at least one component, 

wheiein the adapter translates the output of the at least one component fiom a 
10 fii'st format to a second format; and 

wher ein any other component receiving the output of the at least one component 
receives the translated value.. 

1 1 The method of claiin 2, further compiising; 

(c) coimecthig an adapter' to at least one component, 
15 wherein the adapter applies a function to the output of the at least one component 

to map the oulput fiom a first value to a second value; and 

wherein any other component r eceiving the output of the at least one component 
r eceives the mapped value . 

12 The method of claim 2, fiuthei' compiising: 

20 (c) connectuig an adapter to at least one component^ 

wheiein the adapter tr anslates the input of the at least one component from a fir st 
format to a second format; and 

wher ein the component receives the translated value as input, 
13 . The method of claim 2, further comprising: 
25 (c) cormecting an adapter' to at least one component, 
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whet ein the ad^ter applies a fimctiou to the input of the at least one component 
to map the input fiom a first value to a second value; and 

wheieia the component receives the mapped value as input. 

14. A method in accordance witii any of claims 10 thiough 13, fuither con^iising: 

5 (d) determining whethei an output vaiiable of a fiist component attached to 

an input vaiiable of a second component can supply to the second component a valid 
value to tiie second coniponent; 

(e) if the fiist component cannot supply a valid value, seaicbing in a database 
of adapters foi' an adapter that can be connected to either the first component oi* tiie 

10 second component so as to tiansfoim the value of the output vaiiable of the fiist 
component to a value that is valid for the second con^onent; and 

(f ) if an adaptei' is found, i^etiieving the adapter and connecting it to the fust 
component oi' the second component, 

15 . A method of modeling a biological system, compiising: 

1 5 connecting each of a plurality of instantiated software components to at least one 

othei of the plurality of software components, 

wherein data are diiectly communicated only between softwaie components 
diiectly connected to each othei , 

1 6 . A system foi modeling a biological system, comprising: 

20 means foi' instantiating a pluialily of software components; and 

means foi connecting each of the plmality of software components to at least one 

othei^ of the pluialily of softwaie components, 

wheiein data ai^ diiectly commnnicated only between software components 

directly connected to each othei'. 
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17. A computer-readable mediiam having stoied thei^eon computer-executable 
instmctions foi pezfoiming the steps con5)iismg: 

(a) instantiating a plurality of software components; and 

(b) connecting each of the pluiality of softwaie components to at least one 
5 othei^ of the plmality of software con^onents, 

whetein data are diiectly communicated only between software components 
directly comected to each othea, 

1 8 A system foi modeling a biological system, comprising: 
a processoi; and 
10 a memoiy in commimication with said processoi^ 

wherein said piocessoi causes a pluiality of softwaie components to be 
instantiated in said memoiy; 

wh^ein said processoi causes each of the pluialiiy of sofhvare components to be 
connected to at least one othei of the pluiality of software components; 
15 whetein data are diiectly communicated only between software components 

directly connected to each ofher. 

19 . A method of modeling a biological system, compiising: 

(a) defming and instantiating a software component dynamically; and 

(b) subsequently ledefimng the softwaie component. 

20 20 . The method of claim 19 wheiein the software component is defined by prototype. 
21., A computei-readable medium having stoied theieon a soflwaie component 
l elatiag to a model of a biological system, compiising: 
an attribute relating to at least one state vaiiable; 
an attribute relating to an initialization method; and 
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an attribute relating to a method of detennining the late of chaage ovei time of 
the at least one state variable. 

22, The software component of claim 21, wherein the at least one state variable 
compxises a state variable relating to the concentiation of a chemical,. 
5 23 . A computer -leadable medium having stored thereon an adapter* relating to a 
model of a biological system, comprisiag: 

an attzibute relating to a fust data format; 

an attiibute relating to a second data format; and 

an attribute relating to a function for translating data fiom the first data format to 
10 the second data format, 

24. A computer-readable medium having stored flieieon an adapter relating to a 
model of a biological system, comprising: 

an attiibute relating to an expected output variable of a fir st component; 
an attribute relating to an expected input variable of a second component; and 
15 an attribute relating to a function for mapping the value of the output vaiiable to a 

value that is a valid value for the input variable,. 

25 . A method of evaluating at least one chaiacteristic of a modeled biological system, 
compiising a plur ality of connected sofbvaie components, the method comprising: 

(a) estimating the amount of time necessary to evaluate the at least one 
20 characteristic on a first computer; 

(b) if the estimated time exceeds a predetermined length of time, selecting a 
plmality of connected software components that can be evaluated within tiie 
predetermined length of time on another" coinputer- and that have not previously been 
selected in any iteration of step (b); 
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(c) if any software components have been selected in the cmient iteiation of 
step (b), transmittiag to the ofhei computei- data relating to the attributes of the selected 
softwai^j components, data indicative of the state of fee selected software components, 
data indicative of the late of change of the state of the selected software components, and 

5 a direction to evaluate the selected software components without refeience to the othei 
software components in die modeled biological system; 

(d) r epeating steps (a) thiough (c) with respect to the modeled system without 
any softwaio components pxeviously selected in any iteiation of step (b) until the 
estimated time does not exceed the predeteimlned lengfli of time; 

10 (®) evaluating any softwaxe components not selected ia any iteiation of step 

(c) without lefeience to the othei software components in the modeled biological system; 

(f) iweiving fiom each othei computer to which data was transmitted iq any 
iteiation of step (c), data indicative of a rate of change; and 

(g) evaluating the entiie modeled biological system based on data calculated 
15 in step (e) and received in step (f) , 

26 The method of claim 25, wheiem the data relating to the attiftjutes of the selected 
software components transmitted to Ihe other computer in step (c) are transmitted usitig a 
platform independent distributed component-based aichitectuie.. 

27., The method of claim 25, wherein step (g) comprises (1) dividing the evaluation 
20 of the entire modeled biological system into a plmality of tasks, (2) peiformiug at least 
two of the plurality of tasks in parallel, (3) combining the results of the performed 
parallel tasks, and (4) performing any of the plurality of tasks not performed in parallel, 
28.. The method of claim 27, wherein the tasks performed iq parallel are performed 
on separate processor s.. 
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29. The method of claim 27, wheiein the tasks peifoimed in parallel at© performed 
on separate computeis- 

30, The method of claim 27, wheiein the tasl^ peifoimed in paiallel aie peifoimed 
on separate con^uters liiiked together by the Intemet. 

5 31. The method of claim 27, wherein the tasks peifoimed in paiallel aie peifoimed 
on sepaiate computeis itmning different opeiatiag systems. 

32. A compntei-ieadable medium having stoied thereon computei-executable 
instructions foi peifoiming the steps con^rising: 

(a) estimating the amount of time necessary to evaluate the at least one 
10 chaiacteiistic on a fiist computei; 

(b) if the estimated time exceeds a predetermined length of time, selecting a 
plmality of connected software components that can be evaluated within the 
predeteimined length of time on anothei computer and that have not previously been 
selected in any iteration of step (b); 

15 (c) if any software components have been selected in the cunent iteration of 

step (b), transmitting to the othei computei data relating to the attributes of the selected 
softwaie components, data indicative of the state of the selected software components, 
data indicative of the rate of change of the state of the selected software components, an4 
a diiection to evaluate the selected software components without refeience to the othei- 

20 software components in the modeled biological system; 

(d) 1 epeatiag steps (a) through (c) with respect to the modeled system without 
any softwaie components previously selected in any itemtion of step (b) until the 
estimated time does not exceed the pi edeteimined length of time; 

(e) evaluating any software components not selected in any iteiation of step 
25 (c) without lefeience to the othei" software components in the modeled biological system; 
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leceiving fiom each othet computei' to which data was ttBnsmitted in any 
iteration of step (c), data indicative of a rate of change; and 

(g) evaluating the entire modeled biological system based on data calculated 
in step (e) and leceived in step (f).. 
5 33 A method of evaluating at least one chaiacteiistio of a modeled biological system, 
comprising a phnality of connected software components, the method comprising steps 
foi; 

(a) estimating the amount of time necessaiy to evaluate the at least one 
characteristic on a first computer; 
10 (b) if the estimated time exceeds a predetermined length of time, selecting a 

pluiahty of connected softwaie components that can be evaluated within the 
predetermined length of ^ time on another computer' and that have not prwiously been 
selected in any iteration of step (b); 

(c) if any softwaie components have been selected in the current iteration of 
15 step (b), tiansmitting to flie other computer data relating to the attributes of the selected 
software componentSj data indicative of the state of the selected sofiwaie components, 
data indicative of the rate of change of the state of the selected softwaie components, and 
a direction to evaluate the selected softwaie components without reference to the other 
softwaie components in the modeled biologicai system; 
20 (d) repeating steps (a) through (c) with respect to the modeled system without 

any softwaie components previously selected iu any iteration of step (b) until the 
estimated time does not exceed the piedetermined length of time; 

(e) evaluating any software components not selected in any iteration of step 
(c) without reference to the other software components in the modeled biological system; 
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(f) receiving fiom each other computei to which data was transmitted in any 
iteration of step (o), data indicative of a late of change; and 

(g) evaluating the entiie modeled biological system based on data calculated 
in step (e) and leceived in step (f) 

5 34 A system fox evaluating at least one chaiacteristic of a modeled biological 
system, comprising a plmality of connected sofhvaie components, the system 
comptising: 

means foi estimating the amount of time necessary to evaluate the at least one 
chaiacteiistic on a & st computei; 
10 means foi i-epeatedly selecting a plurality of connected softwaie components that 

have not been pi-eviously selected and that can be evaluated within the predeteimined 
length of time on anothei computei' until all of the unselected components can be 
evaluated on the fir st computer within the pi edeteimined length of time; 

means foi transmitting to each ofhei computei data relating to the attiibutes of 
15 any selected softwait? components, data indicative of the state of the selected softwaie 
components, data indicative of the late of change of the state of the selected softwaie 
components, and a direction to evaluate the selected softwai'e components without 
reference to the other softwai e components in the modeled biological system; 

means foi evaluating any unseleoted software wifliout reference to the other 
20 softwaie components in the modeled biological system; 

means foi leceiving fiom each othei^ computei to which data was tiansmitted, 
data indicative of a late of change; and 

means foi" evaluating the entue modeled biological system based on the 
calculated data and the received data. 
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35, A system foi evaluating at least one chaiacteiistic of a modeled biological 
system, compiising a ptuiality of coimeoted software components^ the system 
comprising: 

a fkst proc^soi; 

5 a memoiy in communication with said fiist processoi; and 

at least one additional processoi ia commimication with said first piocessoi, 
wheiein said fiist processor estimates the amount of time necessary to evaluate 

the at least one chaiacteiistic on a JSist computex; 

wheiein, if ' the estimated time exceeds a predeteimined length of time, said fiist 
10 piocessoi selects a pluiality of connected softwaie components that can be evaluated 

within the piedeteimined length of time on on© of said at least one additional prccessoi 

and that have not previously been selected; 

wheiein, if any software components have been selected, said first processoi 

causes to be tiansmitted to the one of said at least one additional processoi data lelating 
15 to the attiibxrtes of the selected softwaie components, data indicative of the state of the 

selected software components, data indicative of the rate of change of the state of the 

selected softwaie components, and a direction to evaluate the selected softwaie 

components without lefeience to the othei software components in the modeled 

biological system; 

20 wherein said first processor evaluates any unselected software components 

without leference to the othei' software components in the modeled biological system; 

whei'oin said fiist processor receives fiom the one of said at least one additional 
processor, data indicative of a rate of change; and 

whei'ein said first processoi evaluates the entiie modeled biological system based 
25 on the leoeived data and the calculated data. 
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36. A system foi modeling a biological system, compiising 

a plurality of softwai e components, each of which is connected with at least one 
othei software component, 

wherein each of said components directly coromunicates only with parent or child 
5 components. 

37 . The system of claim 36, wher ein said components are distributed components in 
a distributed computing enviionment. 

38. A system for modeling a biological system, comprising 
a user' interface running on a first computei^ 

10 a simulation engine tunning on a second computer , 

wherein said user interface and said simulation engine each comprise at least one 
distributable component. 

39. The system of claim 38, wherein said distributable components are CORBA 
components,. 

15 40,. The system of claim 38, wherein said simulation engine inchides an ODE solver 
component. 

41, The system of claim 38, wherein said simulation engine includes a solver fectory 
component 

42 The system of claim 38, wherein said simulation engine includes a model factory 
20 component. 

43 The system of claim 38, wherein said nsei interface comprises a graphical user 
interface.. 

44., The system of claim 38, wherein said user interface comprises one oi' mor e of the 
following components: an XML parsing component; a mathematical equation 
25 generation component; and a visualization component . 
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45. The system of claim 38» wherein said user- inteiface comprises at least two 
components running on at least two diffeient computers. 

46. The system of claim 38, wheiein said simulation engine comprises at least two 
components running on at least two different computers. 
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